Writhing Photons and Berry Phases in Diffusive Wave Scattering 
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We study theoretically the polarization state of light in multiple scattering media in the limit of 
weak gradients in refractive index. Linearly polarized photons are randomly rotated due to the 
Berry phase associated with the scattering path. For circularly polarized light independent speckle 
patterns are found for the two helical states. The statistics of the geometric phase is related to the 
writhe distribution of semiflexible polymers such as DNA. 



Multiple light scattering techniques jjj allow the exper- 
imentalist to probe deep into opaque, strongly scattering 
samples giving information on both static and dynamic 
correlations. The preferred way of interpreting these ex- 
periments is via a field of scalar photons Q which satisfy, 
in the simplest analytic treatments, a diffusion or trans- 
port equation. The phase of the wave associated with 
each path is just proportional to the path length sampled 
by each photon. In this letter we show the importance of 
geometric phases in the evolution of photons in multiple 
scattering which have as their origin the vector nature 
of light. We give explicit expressions for the probability 
distribution of this extra phase in transmission geometry. 
Rather surprisingly, whilst the basic propagation laws for 
light in inhomogeneous media go back over 60 years to 
the work of Rytov |3],[I| , the full implication of his results 
has not yet been exploited in diffusive wave scattering. 
In this letter, we consider the consequences of Rytov's 
observations on the scattering of circularly and linearly 
polarized light. 

We first summarize our main results before giving a 
more extended discussion: For circularly polarized light 
the photon helicity is a conserved quantum number in 
systems with weak inhomogeneities j^j within the eikonal 
approximation. Each possible, multiple scattered, path 
through a sample is highly tortuous and thus has asso- 
ciated with it a writhe </>j/27r A Berry phase, e^i, 
then adds to the simple geometric optical path tpi = q£i, 
calculated in a scalar theory, where q is a wavenumber, 
li the optical path length and = ±1 is the photon 
helicity. In a transmission geometry this leads to two 
possible speckle patterns for a given arrangement of scat- 
ters as a function of the polarization state of the light. 
For linearly polarized light the situation is more subtle: 
For each individual path the geometric phase rotates the 
plane of polarization . However the final experimen- 
tal measurement is a measure of intensity, summing over 
all possible paths. This sum leads to a final polarization 
state which is in general elliptical rather than linear. We 
shall thus characterize the evolution of a linearly polar- 
ized state as a distribution of polarization states on the 
Poincare sphere. Note that a close analog of Rytov's re- 
sult has been studied in the propagation of polarized light 



along a fiber optic: A beam of light transmitted along a 
fiber Q acquires a non-trivial phase when the tangent 
of the fiber t(s) plotted on a unit sphere encloses a non- 
zero area. This experiment is in many ways the simplest 
realization of the Berry phase and is clearly analogous 
to the propagation of a single ray in a weakly scattering 
medium. 

We note that there is an important qualitative differ- 
ence between the evolution of circular and linear polar- 
ization states in weakly scattering media. For the former 
the helicity can be preserved over arbitrarily large dis- 
tances, whereas we shall conclude that the Berry phase 
associated with linearly polarized light always leads to a 
state of random polarization. We thus suggest that the 
natural setting for experimentally studying polarization 
effects in multiple scattering media is in critical, opales- 
cent samples where the gradual density gradients allow 
an exact formulation of the polarization statistics with- 
out backscattering from interfaces. 

Technically we treat the problem of summing over pho- 
ton paths via a mapping onto a semiflexible polymer, 
which treats multiple weak scattering of photons as an 
angular diffusion process. The Berry phase is calculated 
from the writhe j7j of the photon path using methods 
introduced to study the statistical mechanics of DNA 
and other stiff molecules. Our considerations also link 
up with remarks of Js) who gave a local argument for 
the evolution of the polarization vector, in the context of 
multiple scattering, equivalent to that of Berry's. 

Since the DWS technique is very often used to study 
the properties of colloidal systems we start by explaining 
how some of the above ideas arc applied to such multi- 
center scattering systems, where the mapping onto a tor- 
sionally rigid semiflexible polymer is particularly direct 
and simple to understand. In the theory of DWS from 
colloidal samples one defines two characteristic distances. 
The first, £, is the distance between two collisions be- 
tween a photon and a scattering center. The second, £*, 
measures the distance over which a photon must travel 
in order to forget its initial direction of propagation. In 
a simple analogy with stiff polymers one can consider 
that the length £ corresponds to a monomer size while 
£* is equivalent to the persistence length of the polymer. 
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In strongly scattering media, with high contrast between 
the refractive indexes between inclusions and background 
these two lengths are comparable. However by using par- 
ticles large compared with the wavelength of light and 
with low contrast between the dielectric properties of the 
two media we can easily find samples for which I* /£ is 
of order ten. In recent experiments Q in large droplet 
helium aerosols it is has proven possible to increase £* jl 
to over 200. Motivated by this last experimental system 
we present here arguments as to the polarization statis- 
tics of light scattered in a regime of intermediate sample 
thickness, L, such that £ <C L < I* in transmission geom- 
etry. In such samples photons are scattered many times, 
but are still propagating largely in the forward direction. 
This is exactly the geometry studied experimentally in 

We now derive the mapping which allows us to trans- 
pose results known from the statistical mechanics of stiff 
polymers, including now the torsion as well as bending 
degrees of freedom. Following || we shall model the 
medium as an ensemble of randomly oriented interfaces 
neglecting the spherical structure of the aerosol. This 
was shown to give a good qualitative description of the 
scattering statistics. Consider the collision between a 
photon and a single interface. The photon can either be 
reflected, and thus deviated by a large angle, or refracted 
by a small angle which can be calculated using geometric 
optics. The probability of reflection, for a typical impact 
parameter, is comparable to 1Z — (n— l) 2 /(n+ 1) 2 where 
n is the refractive index of the inclusions relative to that 
of the background medium. If we write e = (n — 1) then 
we see that the intensity of the direct beam should decay 
over a length which is at most £ re f = £/e 2 due to back 
scattering. A background of photons scattered through 
large angles is indeed observed experimentally, in addi- 
tion to the main forward beam |^]. 

Most of the beam is refracted at the interface. For 
typical impact parameters the angle of deviation should 
be comparable to e. Since successive deviations of a pho- 
ton are independent this implies that the direction of 
propagation of the photon diffuses as it propagates into 
the medium with angular diffusion coefficient D ~ e 2 /£. 
Thus photons which are not directly reflected at a sur- 
face also turn over a length comparable to £/e 2 which is 
thus our estimate of £* in this weakly scattering limit. 

Let us now consider the evolution of an incident lin- 
early polarized beam. The reflectivity of an interface is 
a function of the plane of polarization with respect to 
the surface; the transmitted beam has a modified polar- 
ization state. In the hypothetical case of perfect trans- 
mission the plane of polarization of the refracted beam 
evolves by parallel transport djnj . Due to the reflections 
the transmitted amplitudes of the two polarization com- 
ponents (defined relative to the local surface orientation) 
are comparable to 1 — 0(e 2 ). This leads to a rotation 
of plane of polarization of the transmitted light by an 



angle ±0(e 2 ) compared to a parallel transported state. 
In the analogy with a stiff polymer this corresponds to 
excitation of a torsional mode. In stiff polymers one can 
define two independent persistence lengths for bend and 
torsional degrees of freedom, l p and l t . These lengths 
are usually comparable. For the case of multiple light 



scattering we see that £ p 



If 2 whereas t% ~ £/e 4 



The "torsional" degree of freedom for photons is frozen 
out at low dielectric contrast as It ^> l p , so that as e — > 
parallel transport of the transmitted component becomes 
exact, in agreement with Rytov's treatment of Maxwell's 
equations in the eikonal approximation. We have thus 
extended the analogy between stiff polymers Jll[] and dif- 
fusing photon paths by studying the torsional degrees of 
freedom of the polymer, which are entirely analogous to 
the polarization states of photons. 

The tangent to the curve t(s) is a unit vector and can 
be considered as living on a unit sphere. We now use 
the Berry formula linking the area enclosed by the curve 
t(s) on this sphere with the geometric phase to calculate 
the probability distribution of phases associated with a 
samples in the limit I <C L <C I* ■ The orientation of a 
photon diffuses |To| as a function of curvilinear distance, 
so that the probability of orientation obeys the equation 



(i) 



with V 2 the spherical Laplacian and s the optical path 
length. When the sample is thin compared with £* the 
photons do not diffuse very far from their original propa- 
gation direction. We can thus calculate the phase contri- 
bution from the Berry phase by looking at the problem of 
diffusion on a local, planar approximation to the sphere. 
The probability distribution for the area of random loops 
on a plane was solved by Levy [ fH)|Jl2]| with the result for 
the probability distribution for the Berry phase: 
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2DL cosh 2 ((j)/DL) 



(2) 



From this expression we find the mean square phase as 
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We proceed by studying the Jones vectors describing 
the electric field of a coherent light source. Vertically and 
horizontally polarized light is described by the vectors 
j v = (1, 0) and jh — (0, 1) respectively, whereas circu- 
lar light corresponds to the vectors j± = l/y/2 (l,±i). 
We are ultimately interested in the intensities of the var- 
ious polarization states of the transmitted beam which 
are most easily visualized via the Poincare sphere, fig 
(1). The three axes correspond to a series of measure- 
ments h = Io — /go, i2 = I45 - I-A5 and i 3 = I + — J_. 
Here, io,±45,90 are the normalized intensities measured 
with with a linear polarizer inclined at the subscripted 
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FIG. 1. The Poincare sphere describing the polarization 
state of the light. Linearly polarized states correspond to 
points along the equator while the north and south poles, 
labeled j+ and j- correspond to circularly polarized light. 
Other states are elliptically polarized. Multiple scattering of 
the initial state j v describing linearly polarized light leads 
to a circularly symmetric distribution of polarization states 
centered on j v in samples where L/l* < 1 



angle and I± is the intensity measured with circular an- 
alyzers. 

As stated above for circularly polarized light e^, can 
simply be added to the phase ipi so that the transmitted 
Jones vector can be written in the form 

e-^-L^.) (4) 

for the incident state j±. We now sum over all possible 
paths to find the total scattered amplitude: The sum of 
the individual circularly polarized states also gives rise 
to circularly polarized light as the final state. When 
L/l* ~ 2 the geometric phase introduces a relative phase 
2(/)i between the right and left states comparable to n. 
We thus understand that the speckle pattern for the two 
helical states is similar for very thin samples, whereas for 
thicker samples we find two independent intensity distri- 
butions for each circularly polarized state. 

If we illuminate with linearly polarized light j v = 
(1,0), each transmitted photon is described by its writhe 
cf>i, and the total phase ipi so that the transmitted state 
is 

e# VcosjA ^ J 1 \ (5) 

where we have specialized to samples with L/£* < 1. -0 is 
calculated from the statistics of longitudinal fluctuation 
of a stiff polymer by writing the path length as A = 
/ dsy/l + (dr±/ds) 2 . with rj_ the deviation of the path 
from a straight line. By expanding the square root and 



studying the correlation function (A 2 ), |L3| one finds that 
the fluctuations of path length, £i, for crossing a sample 
of thickness L are of order L 2 /£* , when L < £*. Thus ipi 
is very large, so that ipi mod(27r) is very nearly uniformly 
distributed. 

If we sum the amplitude eq. (|B|) over a large number 
of independent paths we see that we find a vector of the 
form 

A relative phase between the two components of the vec- 
tor jf develops because the random sign of <pi . For small, 
fixed (po the vector jf describes a circle of radius 2(po on 
the Poincare sphere as ipi varies between and 2ir. We 
note, in passing, that the full joint distribution function 
of the writhe and path length, V((p, ip), is closely linked 
with the force-torsion response curves measured in DNA 
micromanipulation experiments jnj 

Thus A and A(p are random variables with <p$ ~ 
y/ {(p 2 ) ~ L/£*, eq. ||). tpi is again uniformly dis- 
trubuted. For thin samples this random Jones vector 
is distributed in a circular disk, centered on the initial 
vector j v , fig(l). The combination of a Berry phase com- 
bined with the widely distributed optical path leads to a 
state of elliptical polarization. For thick samples we find 
a (non-uniform) distribution of vectors over the whole 
Poincare sphere leading to loss of memory of the initial 
polarization state in samples thicker than a few times 
£*. Detailed simulations ]T^ , p^ | have been performed 
to study the decay of polarization in colloidal systems 
in transmission geometry with parameters corresponding 
to polystyrene beads in water. It was indeed observed 
numerically that the evolution of the polarization state 
proceeds by the formation of a circular distribution on 
the Poincare sphere. 

The simulations in JlJJljl also found that the polariza- 
tion state of circularly polarized light decays in thick sam- 
ples. The authors defined new lengths £** rc and £*^ ane -, 
while noting that £** rc > £** ane , over which both circular 
and linear polarization decay. We understand the result 
that £** ~ £* even for circularly polarized light as being 
due to back reflection at surfaces. As argued above this 
leads to a decay in intensity (and hclicity flipping) of the 
beam over the length £ re f comparable £* . If we were to 
perform similar experiments on samples with slow varia- 
tions in refractive index, for instance critical, opalescent 
samples we expect that the reflection occurring at inter- 
faces should be suppressed due to critical broadening. In 
this case the helicity of photons is preserved over a length 
which is presumably exponentially large in the eikonal ex- 
pansion parameter 0], Thus £** rc /£* can be arbitrarily 
large. We reiterate that for experiments with linearly po- 
larized light the geometric phase, together with the phase 
ipi due to the varying path lengths implies that memory 



3 



of the linear state must always be lost on a length scale 
comparable to £*. 

To conclude we have reached an understanding as to 
the evolution of the polarization states in multiple scat- 
tering situations via a mapping onto a writhing poly- 
mer. Despite the extreme simplification of our arguments 
the results are in agreement with detailed simulations 
based on a partial wave simulation using Mie diffraction 
theory |I^,|l5|. Our results are valid for coherent light 
sources; the essential step was the combination of the 
amplitudes of the Jones vectors. In scattering with inco- 
herent sources it is rather the Stokes parameters which 
should be combined, presumably leading to multicom- 
ponent transfer theories such as those discussed in |jl6| . 
The contrast between multiple polarized scattering with 
coherent and incoherent sources deserves closer study. 

In our rather crude description of the colloidal regime 
we have neglected diffraction effects, limiting our treat- 
ment to extremely large droplets, however many of the 
qualitative conclusions should hold even for more mod- 
erate particle sizes where diffraction is more important: 
similar conclusions as to the statistics of the photon paths 
were deduced in || using the Born scattering approxi- 
mation. It would be particularly interesting to perform 
simulations in the limit of low optical contrast in order 
to bring out some of the scaling regimes which may exist 
in this limit. 

Finally we have not considered the nature of the non 
Gaussian statistics of the polarization states due to the 
caustics presumably present in weak scattering limits dis- 
cussed in this letter. A final question that we leave 
open is the nature of two time correlation functions of 
the polarization state: As scattering centers move both 
the writhe of the photons paths and the phase associ- 
ated with each path vary. Does this variation contain 
any interesting information on the dynamics of colloidal 
systems? 

We thank Bart van Tiggelen for introducing us to p4| ] 
together C. Derec and P.-E. Wolf for discussions on light 
scattering methods. 
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